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Abstract 

We construct the equation of state (EOS) of nuclear matter using the relativistic mean 
field (RMF) theory in the wide density, temperature range with various proton fractions 
for the use of supernova simulation and the neutron star calculations. We first construct 
the EOS of homogeneous nuclear matter. We use then the Thomas-Fermi approximation 
to describe inhomogeneous matter, where heavy nuclei are formed together with free 
nucleon gas. We discuss the results on free energy, pressure and entropy in the wide range 
of astrophysical interest. As an example, we apply the resulting EOS on the neutron star 
properties by using the Oppenheimer-Volkoff equation. 
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1 Introduction 



Supernova is a spectacular stellar phenomenon happening when massive stars (M > 8Mq) use 
up the nuclear fuel. Due to the gravitational core collapse, supernova explosion may lead to 
the formation of a neutron star or black hole in the central part of the massive star. In order 
to clarify the mechanism of the explosion and the whole phenomena of supernova explosion, 
it is necessary to study them by performing numerical simulations. For this purpose, we need 
the equation of state (EOS) of nuclear matter in the wide density and temperature range with 
various proton fractions, which appear during the collapse and the explosion. 

There have been several EOS worked out so far for supernova study [|I], 0, |^, ^ ^. Most 
of them are constructed based on the non-relativistic Skyrme Hartree-Fock (SHE) framework. 
Although some of these EOS are carefully constructed in the form for the supernova simulations, 
there are many places where the EOS are discontinuous or not available, and the supernova 
simulations stop there due to the lack of the completeness of the EOS table. Sometimes, they 
are not even documented and it is difficult to fix the problem. Lattimer and Swesty worked 
out the EOS in the compressible liquid drop model based on the non-relativistic framework [^]. 
This EOS table is almost the unique one, with which many simulations are made at present. 

Recently, there is a tremendous development in the description of nuclei and nuclear matter 
in terms of the relativistic many body theory ||^, |^, j^. The relativistic Briickner-Hartree-Fock 
(RBHF) theory provided the saturation close to the experimental point in nuclear matter. 
The relativistic effect provides the strong density dependent repulsion over the non-relativistic 
results due to the three-body Z-graph process. The RBHF theory provides also the strong 
spin-orbit splittings, which are the fundamental quantities of the shell model. Based on this 
finding of the RBHF theory, the phenomenological model was developed using the relativistic 
mean field (RMF) theory P, |lO|, |12|. There are many works performed on nuclear ground 
states using the RMF theory. The results compare very well with experimental data on not 
only stable nuclei but also unstable nuclei, which are being obtained using radioactive nuclear 
beams |T^. 

Having success in describing nuclear properties using the RMF theory, it is extremely im- 
portant to work out the EOS for astrophysical purposes. We have to work out consistent 
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calculations for high density nuclear matter, inhomogeneous nuclear matter and finite nuclei 
under the wide range of density, proton fraction and temperature, which appear inside neutron 
stars and supernovae. For this purpose, we study the properties of dense matter with both 
homogeneous and inhomogeneous distribution in the RMF framework adopting the parameter 
set TMl for the RMF lagrangian. The RMF theory with the TMl is known to provide excellent 
properties of the ground states of heavy nuclei including unstable nuclei |jTl], |l4l- The TMl 
was also used for the giant resonances within the RPA formalism and demonstrated to provide 



good descriptions of giant resonances . 

For the description of inhomogeneous nuclear matter, which appears below p ~ Po/3 with 
po being the normal nuclear matter density, we use the Thomas- Fermi approximation |jT6|. We 
assume that the inhomogeneous matter is composed of a lattice of spherical nuclei. At finite 
temperatures, a proton gas, as well as a neutron gas, may be present outside the nuclei. We 
minimize the free energy of the system with respect to the cell volume, proton and neutron 
radii and surface diffusenesses of the nucleus, and densities of the outside neutron and proton 
gases. We compare the calculated free energy with that of homogeneous matter to determine 
the phase of dense matter. This has to be done for each density, temperature and proton 
fraction in the wide range necessary for the astrophysical use. 

This paper is arranged as follows. In Sect. 2, we describe the RMF theory and the TMl 
parameter set. In Sect. 3, we introduce the Thomas-Fermi method to be used with the Wigner- 
Seitz approximation. In Sect. 4, we show the numerical results. In Sect. 5 we show the neutron 
star profiles using the EOS table with the Oppenheimer-Volkoff equation. Sect. 6 is devoted to 
the conclusions and discussions. 



2 Relativistic mean field theory 



We briefiy decribe the relativistic mean field theory here. All the details are to be referred to 
Ref.[P, |ll], |l^ ^]. We start with the lagrangian given by 



RMF 



1 - Ts ' 



(1) 



^^d^ad^a - -m^a - -g2a - -g^a 



where 



W^""" = d^u'^-d'^uj^, (2) 
F^'^ = d^A'^ — 

We adopt the lagrangian with the non-hnear a and uj terms, which was studied by Sugahara 
and Toki [jll|. It was shown that the inclusion of the non-hnear a terms is essential to reproduce 
the properties of nuclei quantitatively and a reasonable value for the incompressibility, while 
the non-linear uo term was added in order to reproduce the density dependence of the vector 
part of the self-energy of the nucleon obtained in the RBHF theory. 

The lagrangian contains the coupling constants, the meson masses, and the self-coupling 
constants as parameters. We adopt the parameter set TMl, which was determined in Ref.[^ as 
the best one to reproduce the properties of finite nuclei in the wide mass range in the periodic 
table including neutron-rich nuclei. They made the least-square fitting to the experimental 
data of the binding energies and the charge radii of proton magic nuclei from Ca to Pb in- 
cluding unstable Pb isotopes. The RMF theory with the TMl parameter set was also shown 
to reproduce satisfactory agreement with experimental data in the studies of the nuclei with 



deformed configuration |Tj] and the giant resonances within the RPA formalism With the 



TMl parameter set, the symmetry energy is 36.9 MeV and the incompressibility is 281 MeV. 



The parameter set TMl can be found in Ref. |TT|, [18 



The derivation of the Euler-Lagrange equations and the following treatment for the EOS of 
homogeneous nuclear matter at finite temperature is the same as the procedure described in 



Ref. . The expressions for physical quantities such as energy density and entropy, which 

are used in the Thomas- Fermi calculations, can be found in [0. Properties of homogenous 
nuclear matter and supernova matter, which contains leptons, at high densities (above po) at 
finite temperature in the RMF theory have been discussed in the same references. 
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3 Thomas- Fermi method for non-uniform matter 



For the density below 10^'^ g/cm^, where heavy nuclei can be formed together with free nucleon 
gas in order to lower the free energy, we perform the Thomas-Fermi calculation based on the 



work done by Oyamatsu |jT6[. In this case, the system is inhomogeneous as modeled by a 
mixture of electrons, muons, free nucleons, and a single species of heavy nuclei. The leptons 
can be treated as uniform non-interacting relativistic particles. They will play no role in the 
minimization of the free energy for the system with fixed proton fraction. So we mainly pay 
attention to baryon contribution in this calculation. 



As in Ref. |TB[ , we assume that each heavy spherical nucleus is considered to be in the center 
of a charge-neutral cell consisting of a vapor of neutrons, protons and leptons. The nuclei form 
the body-centered-cubic (BCC) lattice to minimize the Coulomb lattice energy. 

It is useful to introduce the Wigner-Seitz cell to simplify the energy of the unit cell. The 
Wigner-Seitz cell is a sphere whose volume is the same as the unit cell in the BCC lattice. We 
define the lattice constant a as the cubic root of the cell volume. Then, we have 



Vceii = a^ = Nb/hb, (3) 

where Nb and hb are the baryon number per cell and the average baryon number density, 
respectively. We calculate the Coulomb energy using this Wigner-Seitz approximation and add 



a correction energy for the BCC lattice [0. This correction energy is negligible unless the 
nucleus size is comparable to the cell size. 

We assume the nucleon distribution functions ni{r) for neutrons {i = n) and protons {i = p) 

as 

r r . wn3 

+ n°"*, 0<r <Ri 



rii (r) 



1-1^'- 



nr', R^<r< Rcell 



(4) 



where r represents the distance from the center of the nucleus and Rceii is defined by the 
relation, 

Vceii ^ yi?Lz- (5) 

The density parameters n*" and ?t,°"* are the densities at r = and r > Ri, respectively. The 
parameters Ri and ti determine the boundary and the relative surface thickness of the nucleus. 
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In order to illustrate these parameters, we show the neutron and the proton distributions in 
Fig.l for the case at T = lOMeV, Yp = 0.3, and pb = 10^^'^ g / cm^ . Here we define the baryon 
mass density as pB = Tn^nB with m„ being the atomic mass unit. For the system with fixed 
temperature T, proton fraction Yp, and baryon mass density pb-, there are seven independent 
parameters in the nine variables; a, n^", n°"*, tn-, n^p-, n"^^, Rp, tp. The thermodynamically fa- 
vorable state will be the one that minimizes the free energy density with respect to these seven 
parameters. In zero temperature case, there is no free proton outside nucleus, so the number 
of the independent parameters is reduced to six. 

In this model the free energy density contributed by baryons is given as 

f = {E-TS)/a\ (6) 
where the energy per cell E can be written as 

E = Ebulk + Es + Ec- (7) 
The bulk energy E^uik and the entropy per cell S are calculated by 

Ehuik = / £rmf ( rin (r) , Hp (r) ) rfV, (8) 

J cell 



S = Srmf ( nn (r) , Up (r) ) d r, (9) 

J cell 

where Ermf and srmf are the energy density and the entropy density in the RMF theory as 
functionals of the neutron density Un and the proton density Up. Note that Ermf and sjimf at 
each radius are calculated in the RMF theory for uniform nuclear matter with the corresponding 
densities and n„. As for the surface energy term Eg due to the inhomogeneity of nucleon 
distribution, we take a simple form as, 

Es= f Fo I V (n„ (r) + Up (r) ) ^ d'r. (10) 

J cell 

The paramter Eq = 70 MeV ■ fm^ is determined by doing the Thomas- Fermi calculations of 
finite nuclei so as to reproduce the gross properties of nuclear mass, charge radii and the beta 



stability line as described in the appendix in Ref. |16]. We adopt the same expression for the 



Coulomb energy term Eq as Eq.(3.8) in the same reference |]T6| 
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We perform the minimization of the free energy for each density, temperature and proton 
fraction in the wide range necessary for the astrophysical use. The phase transition from non- 
uniform matter to uniform matter occurs when the free energy density can not be reduced by 
making heavy nuclei. 

4 Equation of state of nuclear matter 

We construct the EOS of nuclear matter in the wide range of the baryon mass density [ps = 
10^ ~ 10^^-^ g/cm^), temperature (T = ~ 50MeV), and proton fraction {Yp = 0.01 ~ 0.5). 
We treat the uniform matter and non-uniform matter consistently using the same RMF theory. 
Hence, all the resulting thermodynamical quantities are consistent and smooth in the whole 
range in this set of the EOS. We note that the EOS of nuclear matter with Yp = 0.0 (i.e. 
neutron matter) for various densities and temperatures in the above range is also calculated 
assuming the uniform distribution in all conditions. 

We discuss first the calculated properties of uniform matter. We show in Fig. 2 the energy per 
baryon of homogeneous nuclear matter at zero temperature as a function of the baryon density, 
which is the basis of this work. We compare the results of the RMF theory (solid curves) 
with those obtained using the parameterized form based on the non-relativistic framework in 
reference ^ (dashed curves). It is shown that the RMF theory provides the strong density 
dependent repulsion over the non-relativistic results. The calculation of free energy of uniform 
matter at finite temperature has been also done in the same manner. Note that the free energy 
at zero temperature is equal to the internal energy. Calculated properties of homogenous 
nuclear matter at finite temperature have been discussed in Ref. The free energy per 

baryon of homogeneous nuclear matter is used as an input in the Thomas- Fermi calculation. 

By making the Thomas-Fermi calculations and comparing the free energies of uniform and 
non-uniform configurations, we determine the most favorable state of nuclear matter at each 
condition of density, proton fraction and temperature. In zero temperature case, the thermo- 
dynamically favorable state for pB < lO^'^g/cm^ with Yp > 0.3 is to have a BCC lattice of 
heavy nuclei without outside neutron gas. When the proton fraction Yp decreases to around 
0.3, the neutrons drip out from the nucleus. In Fig.3 we show the neutron and the proton 
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distributions along the straight hne joining the centers of the nearest nuclei in the BCC lat- 
tice for the case at T = OMeV, Yp = 0.4, and pe = 10^^ 10^^ 10l^ lO^^-^^/cm^ The phase 
transition from non-uniform matter to uniform matter takes place around pb ~ lO^^g/cm^. 
In this paper, we assume only the spherical nuclear shape and neglect non-spherical nuclear 
shapes, which may further reduce the free energy just before the transition to uniform matter. 
However, these shapes would not give siginificant change to the EOS of matter but only make 
the transition smoother. The non-spherical shapes can reduce slightly the surface and Coulomb 
energy compared with the spherical shape, when the size of the nucleus is comparable to the 
cell size as in the case of pb = lO^^g/cm^ in Fig.3. However, in these cases, the surface and 
Coulomb energies are much smaller than the bulk energy so that the energy gain due to the 
shape change is relatively small. We shall calculate the nuclear mass number, A — Z + N, in 
the higher desity part of the inhomogeneous nuclear matter by integrating the proton number, 
Z, and the neutron number, N, up to the radius Rmax, which is the maximum between Rp 
and Rn- Here Rmax is considered as the boundary of the nucleus. In Figures 4a and 4b, we 
display the nuclear mass number A and the proton number Z as a function of the baryon mass 
density Pb with various proton fraction Yp. Both A and Z increase sharply just before the 
phase transition occurs, reflecting that the nucleus becomes larger and the system turns into 
the uniform matter phase. 

In finite temperature case, the second term in Eq.(6) has appreciable contribution to the 
free energy density as the temperature increases. We show the phase diagram of non-uniform 
matter and uniform matter at T = 0, 1, 5 and lOMeV in Fig.5. For extremely low density and 
finite temperature range, the thermodynamically favorable state is uniform nucleon gas. We can 
approximate the uniform nucleon gas as a classical non-interacting gas below pb < 10^^ g/cm^. 
The phase transition from uniform nucleon gas to non-uniform matter takes place at the density 
where the system can lower the free energy by making heavy nuclei. It is seen in Fig.5 that 
the lower phase transition density (from uniform nucleon gas to non-uniform matter) depends 
on the temperature very strongly, while the higher phase transition density (from non-uniform 
matter to uniform matter) is nearly independent of the temperature, which is found almost 
constant at ps ~ IQ^^g / cm^ . As the temperature increases, the area of the non-uniform matter 



8 



phase becomes smaller. The non-uniform matter phase disappears when the temperature is 
higher than ~ 15MeV. We show in Fig. 6 the neutron and the proton distributions for the case 
at pb — 10^^'^ g/cm^, Yp — 0.4, and T — 0,5, 10, IbMeV. From this figure, we can find how the 
properties of the Wigner-Seitz cell change as the temperature increases. The neutron and the 
proton densities outside the nucleus ( and ) become large numbers at high temperature, 
while = at T = OMeV for any ps and Yp. There is a tendency that the nucleon number 
ratio between inside and outside the heavy nucleus decreases as the temperature increases. 
Furthermore, the heavy nuclei are dissociated when the temperature is high enough. 

Fig. 7 shows the free energy per baryon f /ub as a function of the baryon mass density Pb 
with various 1^ at T = 0, 1, 5 and lOMeV. In the uniform nucleon gas phase at low density and 
finite temperature regime, the f /ub are not significantly altered by variations in Yp, because 
the interaction between nucleons is very weak at low density. When the system changes into 
the non-uniform matter phase, the formation of heavy nuclei can reduce f /ub, and brings 
appreciable Yp dependence of f /ub- All curves display a rapid rise in the high density range, 
which is due to the strong density dependent repulsion among nucleons in the RMF theory. We 
show the results of the entropy per baryon as a function of with various 1^ at T = 1, 5 and 
lOMeV in Fig.8. The effect of the non-uniform matter phase on the entropy clearly appear in 
this figure. The drop of the entropy in the non-uniform matter region is due to the formation of 
heavy nuclei. The Yp dependence of the entropy in the non-uniform matter region is opposite 
to the one in the uniform matter regions when the temperature is low. This is because more 
heavy nuclei are formed in non- uniform matter phase at higher Yp, then they can bring more 
rapid drop in the entropy per baryon curves comparing with the case at lower Yp. On the other 
hand, the free nucleon gas outside is a dominant composition of matter at high temperature, 
and the effect of the formation of heavy nuclei is weakened as the temperature increases. 

The pressure can be calculated from f /tib by the following thermodynamical equilibrium 
condition. 



2 dif/ne) 



We show in Fig. 9 the pressure as a function of pb with various at T = 0, 1, 5 and lOMeV. In 
the uniform nucleon gas phase at low density and at finite temperature, the pressure is mainly 
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provided by the thermal pressure, which is given by P = ubT. In this case, the pressure is 
independent of Yp. When the system goes into the non-uniform matter phase, the reduction of 
f/ns due to the formation of heavy nuclei brings a sharp drop in pressure. For high Yp, the 
pressure drops to the negative area. The drop of the pressure in the non-uniform matter phase 
has a strong Yp dependence. In high density region, all curves come back to the positive area, 
and have a rapid rise as the densiy increases. The behavior of the pressure at high density 
is determined predominantly by the contribution of the vector meson, u, and isovector-vector 
meson, p, in the RMF theory [|T^. Because of the non-linear term of the u meson in Eq.(l), the 
repulsive u meson field is suppressed at high density. So it does not increase linearly in density 
as those of another RMF theory without this non-linear u meson term. As a result, the pressure 
is reduced largely at high density as compared with the case in the other RMF theory. We 
stress that this behavior is accord with that in the relativistic Briickner-Hartree-Fock (RBHF) 
theory, on which the RMF theory is based. 



5 Neutron star matter and neutron star profiles 

We construct the EOS of nuclear matter with various proton fractions so that it can be used in 
the simulation of supernova explosion where the /3-equilibrium may not be achieved within its 
timescale. In the study of cold neutron stars, the /3-equilibrium condition can be assumed as for 
static properties. In order to calculate the properties of neutron star matter at zero temperature, 
we have to incorporate the contribution of leptons in the EOS. The proton fraction of neutron 
star matter, which is determined from the /5-equilibrium condition, is shown by the dot-dashed 
curve in Fig. 5. Large proton fractions at high densities beyond l{]^^g/cm? is remarkable in 
the RMF theory [|ll|, [1^, |18[. We display in Fig. 10 the pressure of neutron star matter as a 
function of the mass density. Here the mass density is equal to the total energy density divided 
by c^. Since we have taken account of the non-uniform matter phase, the present EOS named 
NURMF (non-uniform relativistic mean field) is largely different from the EOS named URMF 
(uniform relativistic mean field), which is obtained by assuming uniform matter in the same 
RMF theory for low density range [jT8[. The EOS of BJ obtained by Bethe and Johnson ||T^ is 
also shown for comparison in Fig. 10 by the short-dashed curve. 

10 



Applying the EOS of neutron star matter, we calculate the neutron star profile by solving 
the Oppenheimer-Volkoff equation. The gravitational masses of neutron star as a function of 
the central mass density are displayed in Fig. 11. The neutron star mass is not signficantly 
altered by using the EOS of NURMF comparing with the result using the EOS of URMF. This 
is because the neutron star mass is determinded predominantly by the behavior of the EOS at 
high density. On the other hand, the properties of the EOS at low density are important in the 
description of the neutron star profile in the surface region, where the matter is inhomogeneous. 
In order to illustrate the neutron star profile, we show the baryon mass density distribution of 
the neutron star having a gravitational mass of IAMq in Fig. 12. We display the boundaries 
of the uniform matter and the non-uniform matter with and without free dripping neutrons 
by short-dashed and dot-dashed lines. It is difficult to arrive at the low density region in the 
calculation of the neutron star profile by using the EOS of URMF, and it is also not appropriate 
to describe the neutron star matter by uniform RMF approach as < po/3 where heavy nuclei 
are formed. In Fig. 13, we expand the surface region of the Fig. 12 in order to see the difference 
of the neutron star profiles between NURMF and URMF at the large radius region. It is 
interesting to mention that the neutron star radius increases by having the matter composed of 
nuclei in the most outer part of the neutron star. For the neutron star having a gravitational 
mass of IAMq, the central baryon density is around l.Spo- In this case, the proton fraction in 
the central core of the neutron star is around 0.19 as can be seen by the /3-equilibrium curve 
in Fig. 5. 

6 Conclusions and discussions 

This relativistic EOS of nuclear matter is designed for the use of supernova simulations and 
the neutron star calculations in the wide density and temperature range with various proton 
fractions. We have adopted the RMF theory with the non-linear a and u terms, which was 
demonstrated to be successful in describing the properties of both stable and unstable nuclei as 
well as reproducing the similar self-energies in the RBHF theory. The relativity plays an essen- 
tial role in describing the nuclear saturation and the nuclear structure 0, it also brings some 
distinctive properties in the EOS comparing with the case in the non-relativistic framework. 
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Therefore, it is very interesting and important to study the astrophysical phenomena such as 
supernova explosion, neutron star coohng and neutron star merger using the relativistic EOS 
through extensive comparison with the previous studies with the non-relativistic EOS. 

We have performed the Thomas-Fermi calculation to describe the non-uniform matter which 
appears below pB ~ Po/3. The same RMF theory as the one used at high density has been 
adopted as the input of the Thomas-Fermi numerical code. Hence, all the thermodynamical 
quantities such as free energy, entropy and pressure are found smooth in the resulting EOS. 
The non- uniform matter phase appears in a wide — 1^ region in the phase diagram at low 
temperature, while it completely disappears as temperature increases beyond 15MeV. The 
formation of heavy nuclei in non-uniform matter provides remarkable behavior in the EOS. 

It is also important to pay attention to the experimental data of unstable nuclei, which are 
being obtained using radioactive nuclear beams, so that we can get more information about the 



properties of asymmetric nuclear matter from experiments [^]. It is definitely important and 
meaningful to study unstable nuclei and astrophysical phenomena simultaneously in the same 
framework. 
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Figure captions 

Fig.l. The neutron distribution and the proton distribution in the Wigner-Seitz cell at T = 
lOMeV, Yp = 0.3, and pb = ^^^■^g/cm^ are shown by solid curve and dashed curve, 
respectively. The parameters n*", ?t,°"*, Ri and ti of the density distributions for protons 
and neutrons are defined in Eq.(4), while Rceii represents the radius of spherical Wigner- 
Seitz cell. 

Fig. 2. The energy per baryon as a function of the baryon density for homogeneous nuclear 
matter at zero temperature with the proton fraction 1^ = 0, 0.1, 0.2, 0.3, 0.4, and 0.5. The 
solid curves represent the results of the RMF theory, while the non-relativistic results 
taken from |^ are shown by dashed curves for comparison. 

Fig. 3. The neutron distribution (solid curves) and the proton distribution (dashed curves) 
along the straight lines joining the centers of the nearest nuclei in the BCC lattice for the 
case at T = OMeV, Yp = 0.4, and pe = 10^^ 10^^ W\ W^-'^g/cm^. 

Fig. 4a. The nuclear mass number A as a function of the baryon mass density pb with the 
proton fraction 1^ = 0.1, 0.2, 0.3, 0.4, and 0.5 at zero temperature. 

Fig. 4b. The nuclear proton number Z as a function of the baryon mass density pb with the 
proton fraction 1^ = 0.1, 0.2, 0.3, 0.4, and 0.5 at zero temperature. 

Fig. 5. The phase diagrams of the nuclear matter at T = 0, 1, 5, and lOMeV in pb — Yp 
plane. The regions denoted by Matter correspond to the uniform matter described by 
RMF theory. The regions denoted by Nucleus+Gas and Nucleus correspond to the non- 
uniform matter with and without free nucleon gas, the boundary between them is shown 
by dashed curve at T = OMeV. The regions denoted by Gas correspond to the uniform 
matter close to the classical ideal gas approximation. The /3-equilibrium curve is also 
shown at T = OMeV by dot-dashed curve. 

Fig. 6. The neutron distribution (solid curves) and the proton distribution (dashed curves) 
along the straight lines joining the centers of the nearest nuclei in the BCC lattice for the 
case at pb = lO^''-^g/cm^, Yp = 0.4, and T = 0, 5, 10, 15MeV. 



Fig. 7. The free energy per baryon as a function of the baryon mass density pb with the proton 
fraction Yp = 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5 at T = 0, 1, 5, and lOMeV. 

Fig. 8. The entropy per baryon as a function of the baryon mass density pb with the proton 
fraction Yp = 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5 at T = IMeV, and with Yp = 0.01 and 0.5 at 
T = 5 and lOMeV. 

Fig. 9. The pressure as a function of the baryon mass density pb with the proton fraction 
Yp = 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5 at T = 0, 1, 5, and lOMeV. 

Fig. 10. The pressure of neutron star matter at zero temperature is shown as a function of the 
mass density. The sohd curve is the EOS of NURMF including the non-uniform matter 
phase at low density range, while the dashed curve is the EOS of URMF done by just 
assuming uniform matter phase in the same RMF theory. The EOS of BJ is also shown 
by short-dashed curve for comparison. 

Fig. 11. The gravitational masses of neutron star with the EOS of NURMF and URMF are 
shown by the solid curve as a function of the central mass density. There is no obvious 
difference between the results using the EOS of NURMF and URMF. The neutron star 
masses calculated with the EOS of B J are shown by the short-dashed curve for comparison. 

Fig. 12. The baryon mass density profiles of the neutron star having a gravitational mass of 
1.4M0 using the EOS of NURMF and URMF are shown by sohd curve and dashed curve, 
respectively. We display the boundaries of the uniform matter and the non-uniform 
matter with and without free dripping neutrons by short-dashed and dot-dashed lines. 
There is no obvious difference between the results using the EOS of NURMF and URMF 
in the central core of the neutron star. The difference exists in the surface region as can 
be seen in Fig. 13. 

Fig. 13 The expansion of Fig. 12 in the surface region. The difference between the neutron star 
profiles with the EOS of NURMF and URMF is seen at the large radius region. 



15 




Fig.2 



o.ioh T=OMeV Yp=0.4 P^=^ 



0.05 - 



0.00 
0.10 



-L 



-L 



-L 



_L 



_L 



T=0 MeV Yp=0.4 Pb= 1 0^ ^g/cm^ 



0.05 



0.00 



0.10 - 




T=0 MeV Yp=0.4 Pb= 1 0^ 



0.05 

0.00 
0.10 





-L 




T=0 MeV Yp=0.4 Pb= 1 0^ 



0.05 



0.00 




20 



40 



60 

r (fm) 



80 



100 



120 



Fig.3 




Fig.4a 




6 8 10 12 14 



Fig.4b 



0.5 

.2 0-4 
o 

2 0.3 



§ 0.2 



0.1 - 



0.4 
0.3 
0.2 

1 - T=5 MeV 



0.4 
0.3 
0.2 
0.1 

0.4 
0.3 
0.2 
0.1 
0.0 



Gas 



T=10 MeV 

I L 



Nucleus+Gas 



Matter 



\J 



Gas 



Gas 



- T=1 MeV 



Nucleus 

+ 

Gas 



Matter 



Nucleus+Gas 



V 



P Equilibrium Curve 



Nucleus 



Nucleus+Gas 



T=0 MeV 

I 



Matter 



Matter 



1 



1 2 



1 4 



logio(pB) (9/cm ) 



Fig.5 



0.10 



13.5 , 3 



T=15MeV Y=0.4 Pg=10 g/cm 



c°" 0.05 



0.00 
0.10 



-L 



-L 



-L 



T=10MeV Yp=0.4 Pb=1 



0.05 

0.00 
0.10 




T=5MeV Yp=0.4 Pg=1 



0.05 

0.00 
0.10 



0.05 



0.00 




T=0 MeV Yp=0.4 Pg=1 0^^'^g/cm^ 




1 20 



Fig.6 




Fig.7 




Fig.8 




Fig.9 




Fig. 10 




Fig. 11 




Fig. 13 



